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■ We suggest to perform systematic measurements of the elliptic flow fluctuations which are sensitive 

to the early stage dynamics of heavy-ion collisions at high-energies. Significant flow fluctuations 
are shown to be generated due to the formation of topological clusters and development of the 
filamentation instability. The statistical noise and hydrodynamic fluctuations are also estimated. 
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I. INTRODUCTION 



> 

A high -energy collision of heavy ions is often called the Little Bang because of its similarity to the cosmological Big 
Bang. Both phenomena are violent explosions and both have attracted attention of experimentalists who have gathered 
unprecedented amount of data, limited basically by the data processing technology. The experiments provided a lot 
I of valuable information about the system's evolution. In particular, small variation in the temperature of background 
radiation have revealed mean dipole component, caused by the motion of the Solar System relative to the Big Bang 
""^ ' heat bath. Tiny (~ 10~ 5 ) fluctuations on top of the dipole contribution, which have been recently decomposed into 
angular harmonics with I up to about 2000, show peaks due to frozen sound exited 15 billions of years ago. A study 
1 1 ' of harmonic fluctuations in the Little Bang may possibly reveal something interesting like frozen modulations as well. 

One of the most spectacular experimental results obtained by now in relativistic heavy-ion collisions at RHIC is 
strong elliptic flow quantified by mean value of the second angular harmonics i>2 [1-5]. The phenomenon, which is 
sensitive to the collision early stage [6] when the interaction zone is of the almond shape, is naturally explained within 
a hydrodynamics as a result of large density gradients [7-12]. Since, the hydrodynamic description is applicable for 
a system which is in local thermodynamic equilibrium, the large elliptic flow suggests a surprisingly short, below 1 
fm/c [13], equilibration time which is difficult to reconcile with dynamic calculations, at least those performed within 
the perturbative QCD, see e.g. [14], where the early rapid expansion is closer to free streaming than to hydrodynamic 
evolution. We note here that the hydrodynamic model [15], which assumes the equilibration of only transverse degrees 
of freedom, has appeared rather unsuccessful in describing the experimental data [1-5]. There have also been attempts 
[16,17] to explain the large elliptic flow within models which do not invoke thermodynamic equilibrium. Finally, at 
large transverse momentum px = 2 — 10 GeV the magnitude of i>2 seem to be well described by the surface emission 
model [26]. 

In such a situation, it is certainly desirable to look for experimental observables which can shed more light on the 
system dynamics at a collision early stage. We propose to go beyond measuring the mean elliptic flow magnitude, 
and study its (and higher harmonics, if it ever be possible) fluctuations on the event-by-event basis 1 . To be specific, 
we suggest to measure vi in every collision and then to analyze the variance of v^. The first attempt of such an 
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1 This is similar idea to the above mentioned measurement of the angular fluctuations of T. In cosmology, there is, of course, 
only one event but one can study angular fluctuations in various regions of the sky. 
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measurement has been undertaken in a very recent study by STAR Collaboration [3]. The result is however rather 
inconclusive. Our aim here is to motivate the work in order to improve experimental procedures. 

The elliptic flow fluctuations are shown to be sensitive even to somewhat exotic phenomena which have been argued 
to occur at the collision early stage. We consider here the filamentation instability [18-20] initiated due to the strong 
momentum anisotropy of the parton system, and the generation and subsequent explosions of the topological clusters 
[21]. To detect these dynamical phenomena of interest one needs, however, a reliable estimate of the usual fluctuation. 
Therefore, a magnitude of the statistical noise is discussed in detail. We take into account the interference of various 
Fourier harmonics and finite resolution of the reaction plane reconstruction. We also study the fluctuations caused 
by the impact parameter and particle multiplicity variation. 



II. FORMULATION OF THE PROBLEM 



Two methods have been developed to quantify the ellipticity, as well as higher harmonics, in the azimuthal angle 
distribution. One possibility is to work directly with the correlation functions of 2, 4, 6 or even more particles [22,23]. 
Then, only the relative emission angles of particles matter and there is no need to determine a direction of the impact 
parameter ipR. Using the cumulants, one can partially eliminate detector effects, and consequently even detectors 
with relatively small acceptance can be used. The resulting Fourier coefficients include, however, not only the effects 
associated with the ellipticity of events but arc contaminated by any 2-, 4-, or, in general, n-body correlations caused 
by resonances, jets, quantum statistics etc. However, it is hoped that these nonflow correlations are dominated by 
the two-body effects and thus the genuine four-particle correlations provide rather clean information about the flow. 
Going to the six-particle correlations, the procedure can be further improved. 

In our considerations, we will refer to another method [24,25], which was formulated earlier and is usually called a 
standard one. The method focuses on the angular distributions relative to direction of the impact parameter. The 
experimental procedure splits in two steps which should be as independent from each other as possible. In the first 
step, one uses all available multi-body information about an event in order to determine the impact parameter direction 
ipR. In the second step, one constructs the distribution of the azimuthal angle relative to ipn of 'selected particles' and 
one computes the Fourier coefficients. The sets of particles used at these two steps are different from each other, and 
we will call their numbers as M and N, respectively. In order to reduce non-flow correlations, the particles of both sets 
(subevents) are usually separated by a rapidity gap. They are still correlated by the flow because the direction of the 
impact parameter ipn is a global feature of an event, like the magnitude of the impact parameter itself 2 . In practice, 
it is desired to use the cumulant and standard method simultaneously. Comparing the results, one can eventually 
separate the 'flow' (global correlation between all secondaries) from 'local correlations' involving a few particles only. 
Such a comparison made by STAR [3] has shown that unless one goes to high px or very central collisions, the flow 
dominates and 2-body correlations contribute to V2 obtained from the standard method only at the 10-15 % level. 

Since the one-particle distribution in a single event can be written as 

oo 

P{4>) = ^ [l + 2 v n cos(n{<j> - Vfl))] 0(0) 0(2tt - 0) , (1) 

n=l 

the n— th Fourier amplitude is determined as 



cos(n((f) - ip R )) , 



with 7TT denoting averaging over one-particle distribution in a single event. 

The reaction plane is never reconstructed precisely and the real reaction plane angle ipR deviates from the estimated 
angle ipE- One observes that 



— cos(n(<£ - ip E )) , 



where R n = cos(n(tpn — -0b)) is the reaction plane resolution factor. 



2 In principle, the impact parametr magnitude can be reconstructed from 2-, 4-, 6-, n-particle correlators. However, it 
would be very difficult and presumably rather inaccurate method. Using total multiplicity or forward calorimeter signal does 
the job very well. 
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Let us now think about the ensemble of events with every event representing a single nucleus-nucleus collision. The 
angular harmonics v n are measured for each event. It should be stressed that this is not only the angle ipR (and ipE) 
which varies form event to event but the amplitudes of Fourier harmonics can also vary due to dynamical reasons. 
According to [25], the average over events of the harmonic's amplitude is not defined as 



{v y g UoM<h-i> E )) \ (2) 



but 



;s. Since the procedure of determining the angle -0b is arranged to be 
maximally independent from that of computing of (cos(n{<p — iPe))^, it is expected that the event averaging of R n 

and of cos(n(0 — iPe)) arc independent from each other. If so, 

/ cos(n(0-fe)y \ / 1 \/ , -y \ (cos(n(0-fe))) 

\ — k — / = W/H^ - * b)) / " — m — ' 

where the approximate equality holds if the resolution factor cos(n(V'_R — "0b)) does not much vary from event to event. 
Thus, the definitions (2, 3) are approximately equivalent to each other, provided ipn is reconstructed sufficiently well 
(which requires M » 1). 

While the definition (2) can be uniquely extended to the second moment, there is an ambiguity how to generalize 
the definition (3). Therefore, we define 



cos 



'n(c/>- Vb)) 



{Rn) 

where (• • •) denotes averaging over events. Since the procedure of determining the angle -0b is arranged to be 



and the fluctuations as 

VarK) (O - (v n ) 2 = — ((cos(n(0-^)) 2 ) - (cos(n(0 - fa)) )') , (4) 

where (R n ) 2 enters as a multiplicative factor. However, R n generates event-by-event fluctuations of observed v n . 
To see the effect, let us consider the single-particle azimuthal distribution in a given event of the form (1) with the 
amplitudes v n being exactly the same in all events. Then, (v n ) = v n but 

y , \ _ ( R D ~(Rn) 2 2 



Thus, the fluctuations of R n contribute to Var(w„). We will often use the symbol Sv n = ^/Var(u„). 

In the following sections, we will focus our attention on the second harmonics and consider several sources of the 
v 2 fluctuations. 



III. STATISTICAL NOISE 



We start our discussions of the fluctuations of v 2 with those caused by the finite number N of particles which are 
used at the step 2 of the standard method when the Fourier amplitudes are determined. We assume here that v n do 
not change from event to event. We also assume that the only correlations in the system are those due to the flow. 
Then, the azimuthal distribution of N particles is a product of N single particle distributions. Namely, 

P N (4n,<h, ■■■,4>n)= V n P{<Px) P((f> 2 ) ■ ■ ■ P{<Pn) , 

where Vn is the multiplicity distribution while all distributions P((f>i) are given by Eq. (1). The single particle 
distributions P(4>i) are correlated to each other because of the common angle ipR. 
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In a single event, the ellipticity is found as 

±l£cos(2(^-^)) 



N 

1 1 

V 2 



R 2 N ^ 

where <j>i is the azimuthal angle of i— th particle and N is the event's multiplicity. According to the definition (3) the 
ensemble average of v 2 then equals 

1 / 1 N \ 



i=l 

= VP«T7 #iP(<^»i) / d<hP{<h) ■■ d^ N P{4> N ) V cos(2(& - Vb)) = w 2 • 
iv=i N J o Jo Jo j=1 

We note that the event-by-event averaging of i?2 and of cos(2(0j — "0b)) are assumed to be independent from each 
other. 

The second moment is 

1 / 1 N 2\ 

« = (^ ) (6) 



Tr^^ Vn W d ^ P ^) dfoPifo) ■ ■ ■ / #jvP(«Aiv)(^cos(2(^-^)))' 
" 1 1 ,„,„-,, ,\ / 1 \ /TV - 1' ~ 



Q + (2(^)-l))(^)+« 2 2 (^)(- 



7V 



It has been found in Au-Au collisions at RHIC [3] that (1)4) <C (v 2 ) while (^2) reaches the value of about 0.07 for 
rather peripheral collisions. Taking these numbers into account, we estimate the fluctuations of v 2 as 

where we have also assumed that (N) 1 and that the multiplicity fluctuations are small. 

The second term in r.h.s of Eq. (7) depends on the number of particles M used to determine the impact parameter 
direction. Since M is usually rather large, we assume in accord with [3] that R 2 does not much deviate from unity. 
Then, as argued in [24,25], we have 

(R 2 ) = (cos(2ty* - iM)) « 1 - Mr fa) 2 ) « 1 - ( 8 ) 

where the parameter a depends on the type of weights which are applied. An actual value of a is irrelevant for our 
considerations. Using the arguments which lead us to the result (8), one finds 

{Rf) = (cos 2 (2(^ - Vs))) « 1 - 2<(Vfl - V'e) 2 ) « 1 



(M) • 

Thus, (i?!) — (-R2) 2 ~ (M)~ 2 . Since the number of particles used to determine the reaction plane is larger or at least 
similar to that which is involved in finding v 2 , we conclude that the second term in r.h.s of Eq. (7) can be neglected. 
Thus, we finally estimate the statistical noise as 

Sv 2 = \ . (9) 

(R2)VW) 

As an extra check, we have performed a Monte Carlo simulation of fake events with N particles generated according 
to P{4>) ~ (l + 2u 2 cos(2(</> — iPr))). The obtained variation of v 2 is, of course, in full agreement with the expression 
given above. 

In the subsequent sections, we discuss physical phenomena which lead to the fluctuations of v 2 different than those 
described by Eq. (9), i.e. originating from true event-by-event fluctuations of flow. 
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IV. IMPACT PARAMETER AND MULTIPLICITY FLUCTUATIONS 



As already noted, the observed elliptic flow is naturally described in the hydrodynamic model [7-12]. Therefore, 
are going to discuss here how large are the fluctuations of v 2 within the hydrodynamics. One should distinguish the 
fluctuations due to the varying impact parameter and those due to the thermodynamic fluctuations at fixed collision 
geometry. We start with the former ones. 

As well known, (v 2 ) strongly depends on the collision impact parameter b. In the case of Au-Au collisions at RHIC 
at i/snn = 130 GeV, the dependence has been parameterized [3] as 

(v 2 ) = atb + a 2 b 2 + a 3 , + a 4 b 4 + a 5 b 5 + a 6 b 6 , (10) 

where b is measured in fm and a x = -3.94 • 10~ 4 , a 2 = 2.1 • 10~ 3 , a 3 = -7.06 • 10" 5 , a 4 = -3.2 • 10~ 5 , a 5 = 3.58 • 10~ 6 , 
and cie = —1.17 • 10 . The parameterization assumes that (v 2 ) vanishes for b = and for b > b max = 14.7 fm. 
The v 2 fluctuations due to the varying impact parameter can be estimated by the formula 

Sv 2 = ^Sb. 
db 

The impact parameter fluctuations are, in principle, measurable through the observation of multiplicity of participating 
nucleons N p which in turn is controlled by the collision trigger conditions. Then, 5b can be recalculated into SN P . 
Adopting the linear dependence 

N p = 2z(l--±-), 

y "max ' 

where Z = 79 is the number of protons in a gold nucleus, one gets the formula 

5v 2 = (01 + 2 a 2 b + 3 a\ + 4 a 4 & 3 + 5 a 5 b 4 + 6 a 6 b 5 ) SN p , (11) 

We note that for b w 10 fm where the flow is maximal the v 2 fluctuations due to the impact parameter variation 
vanish because the derivative d(v 2 )/db is then zero. For b = 5 fm where (^2) ~ 0.03, Eq. (11) gives 5v 2 w 8 • 10~ 4 5N P 
which should be compared to the statistical noise (9). For 6 = 5 fm, 5N p = 30 and (N) = 500, the magnitude of the 
v 2 fluctuations caused by the impact parameter variation is approximately equal to that of statistical noise. Thus, not 
only the statistical noise but also the centrality fluctuations must be subtracted from the measured v 2 fluctuations to 
observe the dynamical fluctuations of interest. 

Let us now consider the fluctuations of v 2 due to the variation of thermodynamic parameters. The most important 
are presumably the multiplicity fluctuations. Here, we present some general formulas, considering an example of 
non-statistical fluctuations of multiplicity in the next section. (We follow a similar analysis [21] of the mean px event- 
by-event fluctuations due to radial flow fluctuations.) We assume here that the multiplicity of produced particles is not 
directly used to determine the collision centrality. In such a case the predicted v 2 fluctuations could be significantly 
reduced. 



The fluctuations of v 2 can be estimated as 



which can be rewritten in the form 



where the index P] x (the effective power) is 



- — P, (12) 
(v2)-{N) Ph > (12) 



Ph 



d\n{v 2 ) _ (TV) d(v 2 ) 
dln(N) ~ (v 2 ) d(N) 



with h denoting a hadron species used to determine v 2 . We note that STAR collaboration has already reported the 
data on v 2 for it, K, K si p, A [2]. The first stochastic factor in Eq. (12) is the relative multiplicity fluctuation which 
drives the fluctuations of v 2 , while the second dynamical factor Ph shows how a change in entropy transfers into v 2 . 
Ph is obviously different for various secondary hadron species which can be used to test the idea further. 
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Assuming the poissonian character of multiplicity fluctuations, Eq. (12) can be rewritten as 



Sv 2 = ^LP h . (13) 

The value of the index Ph can be estimated within the hydrodynamics. The calculations presented in [11] for EoS 
LH8 show that changing dN/dy from 200 to 400 leads to the increase of v 2 for pions from 0.028 to about 0.04 in a 
good agreement with NA49 and STAR data, see Fig. 24 in [11]. Reading a logarithmic slope from that figure, we find 
P-n ~ 0.4, which will be used below 3 . 

Comparing Eqs. (9) and (13), one finds that the ratio of the hydrodynamic fluctuations to the statistical noise 
is V2(R 2 }(v 2 )P h which for (R 2 ) = 0.6, (v 2 ) = 0.07 [3] and P h = 0.4 is 0.02. The effect is indeed rather small. 
However, as discussed in Sec. Ill, the magnitude of the statistical noise can be well controlled, and consequently, the 
hydrodynamic fluctuations seem to be detectable. 



V. FLUCTUATIONS INDUCED BY CLUSTER FORMATION 



In general, cluster production induce event-by-event fluctuations of local multiplicity or dN/dy larger than pure 
statistical noise, simply because the number of clusters is smaller than the number of particles. The estimates, which 
we will give, follow ordinary statistical arguments, assuming that the cluster production happens in statistically 
independent way. This is justified by the observation that they all appear at different locations in the transverse 
plane, and also that they depend on specific hadronic and vacuum configurations at the moment of the collision. The 
original presentation of the idea that clusters should lead to observable event-by-event fluctuations of flows has been 
made by one of us at the 2001 CERN workshop [29]. 

Experimentally, existence of strong clustering of produced secondaries in p-p collisions has been known since 2- 
body correlation function was measured at ISR long time ago. To our knowledge, however, precise properties of such 
clusters have been never really explained or well quantified. As a relatively recent example of the cluster study, we 
refer to Fermilab experiment [33] where high multiplicity p-p collisions have been analyzed with the conclusion that 
the average charged multiplicity per cluster is about 4. A single isolated cluster is produced in the so-called Pomeron- 
Pomeron process. Another example is provided by a recent analysis of the old UA8 data [34], showing production of 
clusters of 3-5 GeV mass. It was further shown that such clusters, with mass up to 5 GeV, decay isotropically in their 
rest frame 4 . Understanding of such clusters is very important to clarify a long-standing problem of 'soft Pomeron' 
dynamics. 

As a theoretical motivation, we suggest topological cluster formation in heavy-ion collisions. Inhomogcncous struc- 
ture of the QCD vacuum, with relatively dilute gas of instantons, results in also dilute set of topological clusters 
arising at the collision early stage when the system is promptly excited, from virtual to real classical fields. For 
more discussion of these ideas, specific formulas and original references, see recent paper by one of us [21] and some 
subsequent works [27,28] where the cluster production and decay into gluons and quarks is discussed. For the purpose 
of this paper, it is enough to know that such a cluster, a QCD sphaleron, is like a heavy resonance which is expected 
to decay into about 3 gluons and 6 quarks and antiquarks. In p-p those should hadronizc into specific final states, 
while in heavy ion case these partons are absorbed by the fireball and simply increase the local entropy density. This 
should cause event-by-event fluctuations of radial flow [21] as well as of elliptic flow we discuss in this note. 

Let us now quantify particle number fluctuations caused by the cluster formation. As previously, N denotes particle 
multiplicity in a given px and y window and N c \ is the number of hadrons which can be attributed to the cluster 
decays. With N we denote the hadron multiplicity from all sources different than the topological clusters. Then, the 



3 P. Kolb was kind enough to provide the results for pions and nucleons from his version of hydrodynamics, and the indices 
turned out to be three times smaller, P n w 0.12, P p w 0.13. These number, however, are affected by some artifacts of the 
freeze-out approximations used, especially for low (SPS) energies. As a result, there is a non-monotonous dependence of V2 vs. 
dN/dy with a minimum, reducing the index. We note that a compilation of the AGS-SPS-RHIC data show a monotonous rise, 
as found in [11]. 

4 Unfortunately, the UA2 detector, which was used to collect the data, was just a simple calorimeter, and we do not know 
anything about the structure of these clusters or even mean multiplicities. RHIC detectors and especially STAR can do a lot 
of clarification in p-p mode, provided proper triggers are implemented. 
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relative fluctuation of the particle multiplicity can be written as 



SN = (N ) SN (N cl ) 6N cl 
(N) (N) (N ) + (N) (N cl ) ■ 



Now, we assume that the fluctuations of both the cluster number n c \ and N are poissonian, i.e. 



5n c \ 1 



SN 1 



(«ci> y^M 



(No) ^(Noj ' 



and that every cluster provides k hadrons (N c \ = kn c \). Then, the relative fluctuation of the hadron number is 



where / = (N cl )/(N) denotes the fraction of final state hadrons produced due to the cluster formation. 

Assuming that the cluster production alone is completely responsible for the p-p cross growth with the collision 
energy, one obtains the upper bound 5 on cluster production in Au-Au [35]. Adopting the scaling from p-p to Au- 
Au with the number of hard collisions, it was estimated in [35] that up to roughly 70 clusters per unit rapidity, 
dn c \/dy 70, can be produced around y = in central Au-Au collisions at RHIC. This estimate, in turn, leads to an 
(upper limit) on cluster-related entropy of about half of the total value, and consequently of about half of the total 
multiplicity (/ < 0.5). Keeping in mind that dN/dy i=a 550, one gets k « 4. Inserting these numbers into Eq. (14), we 
find that the formation and subsequent decays of clusters can (maximally) double the multiplicity fluctuations when 
compared to the poissonian fluctuations. 

Using Eq. (12), one immediately translates the multiplicity fluctuations into the fluctuations of V2- Since, the clusters 
can even double SN/(N), the same holds for 5v2/(v2)- Once we have concluded Sec. IV that the hydrodynamic 
fluctuations seem to be measurable, we claim here that there is a chance to detect the fluctuation growth due to 
the cluster formation. However, it should stressed that the production and subsequent decay of the clusters can be 
observed directly studying the multiplicity fluctuations. 



When the momentum distribution of partons is strongly elongated in one direction, say along the z i.e. beam 
axis, the neutral system has a tendency to split into the filaments along z with the current flowing in the opposite 
directions in neighboring filaments. The reason is as follows: once the currents in the system occur, they generate 
the (chromo-)magnetic field, oscillating in the direction perpendicular to the beam axis, and the Lorentz force acts 
back on the charges which form the currents. It appears that the currents get focused and the current magnitude 
grows. This is the filamentation instability [30] which have been studied in the context of ultrarelativistic heavy-ion 
collisions in [18-20]. 

The breakdown of the azimuthal symmetry of the system due to the instability development gives a chance to 
observe it experimentally. It has been argued [18-20] that the instability growth leads to the energy transport along 
the wave vector which coincides with the Poynting vector of the generated chromodynamic field. Consequently, one 
expects significant a variation of the transverse energy as a function of the azimuthal angle. 

Here we point another, presumably more realistic, possibility to detect the color filamentation. When the instability 
grows the trajectories of charge particles are focused in the centers of the filaments. Therefore, according to the 
Liouville theorem the distribution of the momentum perpendicular to the filaments, say along the x— axis, has to 
expand to conserve the phase space volume. The quantum mechanical counterpart of the argument relies on the 
uncertainty relation: once the particles are localized within the filaments their transverse momentum has to widen to 
the inverse filament thickness multiplied by Ti. Below, we quantify this quantum mechanical reasoning. 

It should be clearly stated that the collective motion caused by the instability development is not correlated with 
the reaction plane and it has nothing to do with the hydrodynamic flow. As such it would be called a 'non-flow' effect. 



5 This is an upper bound because nuclear modification of the structure functions is ignored. The realistic number is presumably 
factor 2 or so smaller. 




(14) 



VI. FLUCTUATIONS INDUCED BY FILAMENTATION INSTABILITY 
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However, the filamentation generates a finite value of v 2 and it contributes to its fluctuations. So, let us consider the 
phenomenon in more detail. 

Let the wave vector of the filamentation mode (k) be oriented along the x— axis. Then, the single particle wave 
function describing the transverse degrees of freedom is the form 



ip(x, y) ~ cxp 



y 



AR 2 



cos{kx + a) , 



(15) 



where R is the system transverse radius. To simplify further analysis we put the phase a equal to zero. Then, as we 
will see, odd harmonics of the azimuthal distribution vanish due to the mirror symmetry of the wave function (15). 
Performing the Fourier transform of (15), one gets the momentum distribution as 



P(p x ,Py) = \4>(Px,Py)\ 2 ~ exp[ - 2R 2 (p 2 x +p 2 y )] [cxp(4RVfc) + cxp( - AR 2 p x k) + 2 . 
Since p x = pcoscf) and p y = psin<^>, the distribution (16) provides the azimuthal distribution of the form 



/>oo />oo 

Pifi) = / dpp P(pcoscj),psm(j)) <~ / dpp< 
Jo Jo 



-2R 2 p 2 



^4i^ pkcostp _j_ ^—47? pk cos0 _|_ 2 



(16) 



(17) 



Neglecting the jacobian p in Eq. (17), the integral over momentum can be performed analytically and the result reads 



P{4>) 



2R k cos 2 



+ 1 



(18) 



The distribution (18) gets a particularly simple form when the filament thickness is much smaller than the system 
size. Then, Rk^> 1 and 



P{<f>) = 



l 



2(71-1) 



1 - 5{<j> - tt/2) - 5{<j> - 3tt/2) 



(19) 



Using the distribution (19), one finds 



vi = 



1 



(20) 



The value of v 2 is rather large. A realistic value of v 2 is presumably significantly smaller because not all particles 
produced in a given event would participate in a collective motion caused the instability development. The collective 
motion should be also convoluted with the thermal one. Thus, the effect of filamentation must be diluted. An 
appearance of the instability in the system is not a deterministic but a random process. Therefore, we expect that 
there are collisions with and without the instability. Consequently, v 2 varies between zero and maximal value (20). 
Thus, we expect large fluctuations of v 2 . 

It should be also stressed here that there are specific distinctive features of the collective flow and the flow fluctua- 
tions due to the filamentation. First of all we note that in contrast to hydrodynamically generated v 2 , the flow caused 
by the instability development does not vanish at zero impact parameter. Thus, one should look for filamentation in 
maximally central collisions. It is also expected that particles with small pr are particularly sensitive the collective 
motion of interest. 



VII. FINAL REMARKS 



The aim of this paper is to advocate usefulness of the flow fluctuation analysis in revealing of the early stage 
dynamics of heavy-ion collisions. We have shown that even rather exotic phenomena can be studied in this way. 
Those presented should be, of course, treated only as examples motivating the measurements, and simple order-of- 
magnitude estimates. 

Since v 2 is experimentally determined on the event-by-event basis anyway, the proposed fluctuation measurement 
presumably does not require much additional efforts. However, an accuracy of the measurements should be improved. 
Studying the flow fluctuations as a function of particle multiplicity one can check whether 5v 2 scales like 1 / \J (N), 
which is a characteristic feature of statistical noise. If not we deal with nontrivial dynamical fluctuations. However, 
before such a conclusion is achieved one has to properly subtract the fluctuations due to the impact parameter 
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variation. We note that these fluctuations can be constrained by the collision trigger condition and that the fluctuations 
due to impact parameter vanish around the maximal flow where d{v2) /db = 0. 

To disentangle various fluctuation sources, the data should be analyzed in a broad interval of impact parameters 
and varying acceptance windows. Fortunately, the mechanisms of interest contribute differently to Sv 2 . In particular, 
the cluster effect is expected to be the largest for most peripheral events, while that of the filamentation for the most 
central ones. It would be also desirable to study elliptic flow fluctuations simultaneously with fluctuations of particle 
multiplicity and other collision characteristics. 

At the end, we mention one more supplementary method [31] to study the azimuthal fluctuations which seem to be 
particularly useful to detect non-flow correlations. The method, which uses the so-called <i>— measure of fluctuations 
[32], does not require the reaction plane reconstruction and can be rather easily applied to experimental data. The 
measure is sensitive to various sources of dynamical correlations and the integrated information provided by $ can 
be combined with that offered by the Fourier analysis [24,25]. Since all Fourier harmonics contribute to $ one can 
check whether the measured harmonics saturate the observed value of <F 
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